How do Quasicrystals Grow? 
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We investigate the physical mechanism underlying the formation of the quasicrystal phase from a 
supercooled liquid. Using molecular simulations, we show that quasicrystal growth is controlled by 
a tendency for atoms to retain their local liquid configuration when attaching to the growing solid 
nucleus. In the system under investigation, which forms a dodecagonal quasicrystal, we show that 
this process occurs through the assimilation of stable icosahedral clusters in the supercooled liquid 
by the growing quasicrystal, which hinders periodic crystal formation in favor of aperiodic crystal 
formation. 



PACS numbers: 61.44.Br 61.20.Ja 

Quasicrystals [U are a unique class of ordered solids 
that display long-range aperiodicity, which distinguishes 
them from ordinary crystals. It is not known what "spe- 
cial" qualities systems must possess in order to form qua- 
sicrystals versus crystals. In a basic sense quasicrystals, 
like crystals, form via nucleation and growth 1 , whereby 
a microscopic 'nucleus' of the solid phase spontaneously 
arises in the supercooled liquid and spreads outward, con- 
verting the system from liquid to solid [2 . One of the 
fundamental problems in quasicrystal physics is to rec- 
oncile the typical picture of nucleation and growth with 
the formation of an aperiodic solid. Quasicrystals can- 
not grow like crystals, where the nucleus surface acts as 
a template for copying a unit cell via local interactions. 
Rather, quasicrystals, require specialized "growth rules" 
that dictate their formation|3]. 

Quasicrystal growth rules fall into two categories: 
energy-driven quasiperiodic tiling models |U [S] and 
entropy-driven random tiling models 017]. While energy- 
driven models rely on "matching rules" to dictate how 
atomic clusters or tiles attach to the nucleus, entropic 
models allow tiles to attach randomly to the nucleus with 
some probability. Although these models provide impor- 
tant insight into the rules that systems must follow to 
grow quasicrystals, they do not provide a detailed physi- 
cal grounds to justify why such rules exist. The phys- 
ical driving force underlying quasicrystal growth, and 
whether it is based on local interactions or long-range 
correlations, is not well understood. 

In this Letter, we elucidate the atomistic mechanism 
of quasicrystal growth by studying the formation of a 
metastable dodecagonal quasicrystal from a model su- 
percooled liquid. We show that quasicrystal growth is 
facilitated by structurally persistent atoms that become 
kinetically trapped in the supercooled liquid. As the 
solid nucleus grows, it incorporates these atoms in a way 
that minimizes expensive rearrangements and eases prop- 
agation allowing the quasicrystal to form instead of the 
stable crystalline phase. In the system under investiga- 
tion, we find that structurally persistent atoms are highly 
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FIG. 1: Dodecagonal quasicrystal and approximants. (a) The 
inherent structure of a 17,576 atom dodecagonal quasicrystal 
formed by the Dzugutov system using molecular dynamics at 
T*=0.42 and p*=0.85. The image on the right shows the 
aperiodic tiles formed by connecting the centers of the do- 
decagonal rings of atoms, (b) Unit cells of various quasicrystal 
approximants. 



icosahedral, indicating that they were in icosahedral clus- 
ters prior to nucleation. Our results provide a significant 
step forward in understanding quasicrystal formation and 
are consistent with the class of systems that are known 
to form quasicrystals experimentally. 

To obtain these results, we utilize both conven- 
tional Monte-Carlo (MC) simulations and umbrella 
sampling [8l. First, we use isothermal (NVT) MC to ob- 
serve the growth of the quasicrystal from a static seed 
nucleus. We then use isothermal-isobaric (NPT) MC to 
observe the growth of large quasicrystal nuclei, gener- 
ated via umbrella sampling. Finally, we use umbrella 
sampling to generate many configurations containing nu- 
clei to study the relationship between quasicrystal nuclei 
and icosahedral clusters. All simulations consist of 3375 
atoms with pair interactions modeled via the Dzugutov 
potential [5^ . Upon shifting the energy scale the Dzugutov 
potential tracks the 12-6 Lennard- Jones potential up to 
a distance at which an additional repulsive term domi- 
nates, suppressing the formation of BCC, FCC, and HCP 
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FIG. 2: Dependence of quasicrystal 
tiling arrangement on liquid struc- 
ture. The images show characteris- 
tic results from MC runs with the 
same quasicrystalline seed (blue) but 
with a different random number se- 
quence. At certain points in space, 
highlighted in yellow in the overlay, 
the tiling arrangements differ. 



crystals and favoring polytetrahedral local ordering. This 
causes the system to form a dodecagonal quasicrystal 
from the melt at certain state points [101 E] (see Fig. 
la). Although the quasicrystal is physically stable over 
the timescale of a typical simulation, it is thermodynam- 
ically metastable with respect to the cr-phase periodic 
approximant |10j (see Fig. lb). Here, we run simulations 
at temperature T*=0.55, pressure P*=3.5 and density 
/3*=0.85, which is slightly below the degree of undercool- 
ing at which the system spontaneously forms a quasicrys- 
tal. At this state point it is improbable to form multiple 
nuclei simultaneously. 

We begin by considering the growth of the solid phase 
from a small static seed nucleus in the form of a periodic 
approximant [12j that is inserted into the MC simulation 
cell (see Fig. 2). Periodic approximants are crystals with 
identical local ordering to the quasicrystals; therefore, for 
small nuclei, quasicrystals and approximants are identi- 
cal and the difference in long-range ordering results from 
a different growth mechanism. Constraining the seed al- 
lows us to study the effect of seed structure on the qua- 
sicrystal growth. We randomize the atoms in our isother- 
mal MC simulations at high temperature before quench- 
ing to T*=0.55, at which point atoms begin to attach 
to the seed, causing rapid solidification. We observe that 
the system consistently forms a dodecagonal quasicrystal 
for all seed sizes, positions, and approximant structures, 
indicating that the system does not copy the seed, but 
rather incorporates atoms into the solid via a different 
paradigm. 

Energy-driven quasicrystal growth models suggest 
that atomic attachment to the nucleus is deterministic, 
whereas entropy-driven models suggest that it is stochas- 
tic. We test the applicability of these models for our sys- 
tem by modifying the sequence of random numbers gen- 
erated during the simulation, holding all else constant. 
As depicted in Fig. 2, for the same seed nucleus (blue), 
we consistently obtain distinguishable quasicrystal tiling 
arrangements, indicating that quasicrystal growth in this 
system is stochastic. We note that most of the tiling 
discrepancies (yellow) represent phasons[13], local tiling 
arrangements that can be converted into one another 
via so-called phason flips [71 |T3]. Because the energetic 
difference between phasons is small, phasons represent 
equivalent pathways by which the quasicrystal can grow. 



thereby providing additional degrees of freedom for the 
growth of quasicrystals versus approximants. 

Although the random tiling model appears to provide a 
suitable description of quasicrystal growth in our system, 
we note that discrepancies in the tiles appear to lag be- 
hind changes to the random number sequence. Thus, the 
process is not truly random, but dependent on some sys- 
tem variable that mimics randomness but reacts slowly 
to random number perturbations. The spatial arrange- 
ment of atoms in the liquid phase is a logical candidate 
for such a variable, since it displays both properties. We 
test this directly by detecting structural correlations be- 
tween atoms in the liquid and the quasicrystal tiles that 
they subsequently form. 

To detect such correlations, we introduce two new met- 
rics: an order parameter to differentiate between liquid 
and solid atoms on a per-atom basis, and a configura- 
tional persistence function to measure how closely atoms 
retain their configurations as a function of time. Our 
order parameter is a modification of the qe{i) ■ qeU) 
scheme of reference |14j. Here, q6{i) ■ qeU) measures 
whether two atoms i and j have correlated first neigh- 
bor shells in terms of shape and orientation, with 1 (0) 
representing perfectly correlated (uncorrelated) configu- 
rations. Atoms form a solid-like connection if qe{i) ■ qeU) 
exceeds a certain value, and atoms that form many solid- 
like connections with neighboring atoms are defined as 
being solid-like, reflecting the fact that in simple crys- 
tals all atoms have identical coordination shells. This 
scheme must be modified for quasicrystals and approxi- 
mants, since neighboring atoms have different coordina- 
tion shells corresponding to Frank-Kasper polyhedra[15] . 
For dodecagonal quasicrystals, we increase the range of 
the neighbor cutoff to r^ut — 2.31ct, corresponding to the 
first ^^2.5 neighbor shells. Also, we modify the harmonic 
frequency I to from 6 to 12, to arrive at a new order 
parameter, <7i2(i) • 9i2(j)- Pairs of atoms form a solid- 
like connection if qi2{i) ■ qi2{j) > 0.45. atoms with 50% 
solid-like connections are solid-like, otherwise they are 
liquid-like. These cutoffs are chosen so as to maximize 
the distinction between liquid and quasicrystal. As qua- 
sicrystals and approximants have identical local ordering, 
'Zi2(i) • qi2ij) niakes no distinction between them. 

We use a similar approach to quantify the configura- 
tional persistence of atoms as a function of time. We de- 
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fine an autocorrelation function qQ{t) = qQ{i; to) ■ qQ{i; t), 
where qQ{v, to) ■ qe{i; t) is equivalent to qeii) • 96 (j), except 
that we apply the dot product to the same atom i at 
two difi^erent times instead of to two different atoms at 
the same time. We define Vcut — 1-65 to include the first 
neighbor shell in our analysis. We use I — 6 since this 
yields the strongest signal, although we obtain qualita- 
tively equivalent results for other choices of I. We nor- 
malize qQ{t) on the interval [0, 1], where qe{t) = 1(0) for 
perfectly correlated (uncorrelated) configurations. 

For reference, we consider the ensemble average {qe{t)) 
versus t for atoms in the bulk liquid and the bulk qua- 
sicrystal, shown by the bottommost and topmost curves 
in Fig. 3a, respectively. Both systems sample the 
isothermal-isobaric ensemble at T*=0.55 and pressure 
P*=3.50 corresponding to (p*) = 0.85. We see that 
atomic configurations in the liquid become completely 
decorrelated from their initial state after ~35,000 MC 
cycles when {qeit)) 0, whereas configurations in the 
quasicrystal are persistent over long timescales. The ini- 
tial drop in {qeit)) from 1.0 to 0.7 is due to thermal 
vibrations. 

We detect structural correlations between liquid atoms 
and quisicrystal tiles in many independent nucleation 
events whereby the system grows a quasicrystal. Pre- 
viously, we used a seed to initiate nucleation; however, 
here we use umbrella sampling to generate many config- 
urations with growing nuclei. This allows us to observe 
nucleation at the same moderate degree of undercooling, 
but also ensures that the configurations that we generate 
are physical. Our sampling scheme is similar to references 
^JM [T5] , where our NPT MC runs are biased according to 
the harmonic weight function w = {N — No)) ■ Here, 
k ~ 0.075, N is the number of atoms comprising the nu- 
cleus (measured by 912 (*) • 912 (i)), and A'o is specified 
such that nucleus sizes near No are sampled selectively. 
We slowly increase the bias from A^o = 10,20, ...90 so 
that nuclei reach N = 80 — 100. We then use these mi- 
crostates as starting points for unbiased NPT MC runs. 
We observe that nuclei with N > 75 atoms tend to grow, 
although many factors other than size (e.g., shape, struc- 
ture, etc.) may affect nucleus stability as well|17|. 

We run MC simulations of growing nuclei for 75,000 
MC cycles, the time it takes for nuclei to grow from ~100 
atoms to ~500 atoms. At this point, boundary efi^ects 
due to the finite size of the simulation cell begin to ap- 
pear. The middle curve in Fig. 3a shows {qeit)) vs. t for 
atoms that attach to the growing quasicrystal nucleus at 
to = 0, such that for t < 0, the atoms are in the liquid, 
and for t > 0, the atoms are in the solid nucleus. We 
observe that {qe{t)) exhibits relatively high values for all 
t < 0, indicating that atoms in the quasicrystal are cor- 
related with their former liquid configurations over long 
timescales. 

This interesting behavior may result from structural 
correlations between the quasicrystal and the liquid, or 
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FIG. 3: Structural correlations, (a) Average value of qeit) 
versus t (MC steps). From top to bottom: atoms in the do- 
decagonal quasicrystal, atoms in the non-equilibrium nucle- 
ating system that attach to the nucleus at t = 0, atoms in the 
liquid. For all runs, the reference time to = 0. (b) Probability 
distribution of qeit) for three times t < for atoms that at- 
tach to the nucleus at t = (dashed line). For comparison, we 
show typical distributions where atoms relax normally (solid 
line). The plots are colored based on their corresponding col- 
ored box in (a). 



simply slow relaxation of liquid atoms in the vicinity of 
the nucleus, or both. To ascertain whether quasicrystal- 
liquid correlations play a significant role, we consider the 
probability distribution P {qQ{t)) (see Fig. 3b) for both 
nucleating atoms (dashed line) and atoms for a reference 
state at which the probability distribution exhibits typ- 
ical behavior (solid line, taken here to be the liquid at 
T*=1.0, (p*) =0.85). The distributions are shown at 
time delays for which they are identical for low qeit). We 
observe that atoms attaching to the nucleus a,t t — 
exhibit a noticeably higher proportion of strongly corre- 
lated atoms than the reference state, indicating that the 
quasicrsystal tends to retain the liquid configurations of 
a particular subset of atoms. Thus, we conclude that 
quasicrystal-liquid structural correlations play a signifi- 
cant role in producing the behavior observed in Fig. 3a. 
However, we note that slow relaxation around the nucleus 
is also a factor, since upon truncating the correlated tails 
of the distributions, we still obtain higher {qe(t)) than 
for the liquid at all t. 

It is probably no coincidence that the quasicrys- 
tal tends to grow in a way that allows certain parti- 
cles to retain their liquid configuration. As quasicrys- 
tals typically form in rapidly quenched metallic alloys 
with strong icosahedral ordering [TS] (similar to metallic 
glasses) , atoms in quasicrystal-forming liquids tend to fall 
into local energy minima that resist rearrangement ^19j . 
Such systems should favor quasicrystal formation since 
quasicrystals are seemingly able to reach a "structural 
compromise" with the liquid in order to grow more 
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FIG. 4: Icosahedral environment, (a) Simulation snapshot 
showing a quasicrystal nucleus (red) together with icosahe- 
dral clusters (yellow for icosahedral centers, white for surface 
atoms) in the liquid (blue), (b) The average probability of 
observing a atom at the center of an icosahedron versus r, the 
distance from the nucleus surface. 

freely. In the Dzugutov liquid, atoms organize into lo- 
cal energy-minimizing structures known as "icosahedral 
clusters [5nj HI]," which exhibit lower mobility than the 
bulk [22] • Structural analysis [TS] of the most persistent 
atoms in Fig. 3b (the correlated tails of the distribu- 
tions) shows that they consist of 32% icosahedra, 25% 
Z13, 31% Z14, and 10% Z15 configurations, consistent 
with the structure of icosahedral clusters. The most per- 
sistent atoms exhibit stronger icosahedral character than 
either the liquid (6% icosahedra), or the dodecagonal 
quasicrystal (27% icosahedra). We note that although 
the icosahedral glass formed by the Dzugutov system has 
vibrational modes similar to the (T-phase|53] (25% icosa- 
hedra and 75% Z14), the most persistent atoms do not 
exhibit cr-like character. 

It is curious that our results implicate icosahedral clus- 
ters as the driving force underlying quasicrsytal growth, 
since icoshedral clusters are relatively scarce at the state 
point under investigation |20j . We can obtain a more intu- 
itive picture of the role of icosahedral clusters by consid- 
ering their spatial arrangement in relation to the grow- 
ing quasicrystal nucleus. We generate a large number 
of nuclei using the umbrella sampling scheme outlined 
above. However, here, following reference |16j . we allow 
configuration swapping between simulations via parallel 
tempering. In all, we run 10 simultaneous MC simula- 
tions for 3.5 million timsteps, where each simulation has 
a unique biasing potential minimum iVg = 10, 20, . . . 100. 
for a given simulation. We save configurations every 100 
MC steps, giving us 35,000 total microstates containing 
nuclei of sizes A^=10-110 for analysis. We identify icosa- 
hedral clusters in our microstates using the method of 
reference [2?], an extension of the method of reference 
i25j. 

As depicted in Fig. 4a, we observe a high percentage 
of icoshaedral clusters (yellow, white) in the region sur- 
rounding the nucleus (red). We quantify this by calcu- 
lating Pico{r), the average probability of observing an 
atom at the center of an icosahedron a distance r away 
from the nucleus surface (see Fig. 4b). For nuclei of all 



sizes, we observe that Picoir) starts with a high value of 
0.15 near the nucleus surface and decreases to the Uquid 
value over a range of about three particle diameters. This 
result can be interpreted in one of two ways: either the 
presence of the nucleus encourages the formation of icosa- 
hedral clusters in the surrounding liquid, or vice versa. 
In either case, the growing nucleus must attach atoms 
from low-energy icosahedral clusters to its surface, which 
facilitates quasicrystal growth. 

Our collective results imply that the driving force for 
quasicrystal formation vs. crystal formation is the ability 
for the solid phase to propagate more freely by growing 
aperiodically vs. periodically. In the system under inves- 
tigation, this process is driven by the tendency for atoms 
to become kinetically trapped in icosahedral clusters in 
the supercooled liquid. However, we note that this can- 
not be the case in general since many quasicrystals do 
not contain icosahedral clusters. Despite this, we spec- 
ulate that the basic mechanism at hand — the tendency 
for low-energy atoms in the liquid to retain their configu- 
ration when incorporated into the growing solid nucleus 
— may be widely applicable to quasicrystal-forming sys- 
tems. 
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